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We discuss the relationship between the cycle probabilities in the path-integral representation of 
the ideal Bose gas, oflt-diagonal long-range order, and Bose-Einstein condensation. Starting from the 
Landsberg recursion relation for the canonic partition function, we use elementary considerations 
to show that in a box of size the sum of the cycle probabilities of length k ^ equals the 
off-diagonal long-range order parameter in the thermodynamic limit. For arbitrary systems of ideal 
bosons, the integer derivative of the cycle probabilities is related to the probability of condensing k 
bosons. We use this relation to derive the precise form of the in the thermodynamic limit. We 
also determine the function tt^ for arbitrary systems. Furthermore we use the cycle probabilities to 
compute the probability distribution of the maximum-length cycles both at T = 0, where the ideal 
Bose gas reduces to the study of random permutations, and at finite temperature. We close with 
comments on the cycle probabilities in interacting Bose gases. 



I. INTRODUCTION 

The canonic partition function of N ideal bosons in a 
system with energy levels eo = < ei < . . . at inverse 
temperature /3 = l/(fcBT) is given by 



Zn = 



N 



HE 

. a n„=0. 



exp 



,N- 



(1) 



For each state cr, the allowed occupation numbers go 
from to N, and the Kronecker (5-function enforces the 
total number of particles to be N. 

The Feynman path-integral [1] represents Zn as a trace 
over the diagonal density matrix p^°^, 

Zn = J dxi .. . dxNP^°''{{xi, . . . , xn}, {xi, . . . , xn}, /3) , 

which is given in terms of the iV! permutations P of the 
distinguishable-particle density matrix p 



xn},Wi, - ■ ■ ,x'n}^P) 
^ E ^('t^i' ■ ■ • ' J'' ^^'Pi ' ■ • ■ ' ^Piv }' /3) • (2) 



iV! 



Thus the partition function is a sum of A^! permutation- 
dependent terms, Zn = "^ZpZp, and it is possible to 
define the probability of a permutation P as 



Zn 



(3) 
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Permutations can be broken up into cycles, and one may 
also define cycle probabilities: 



1 

Zn 



= ^ > Z, 



Pe{N.k} 



(4) 



where {TV, fc} denotes the set of all permutations where 
the particle iV belongs to a cycle of length k. The choice 
of the particle iV in eqn (4) is arbitrary. 

In this paper, we are concerned with the character- 
ization of the cycle probabilities tt^ in the ideal Bose 
gas and with their relation to quantities characteriz- 
ing Bose-Einstein condensation, namely the off-diagonal 
long range order (in Section III, see also [2, 3]) and the 
condensate fraction (in Section IV). In fact, for arbitrary 
finite systems of ideal bosons, the discrete derivative of 
the function tt^ yields the probabilities for condensing 
k bosons [4]. This relation between the cycle statistics 
and the condensate fraction has not been scrutinized be- 
fore in detail. It allows to determine the tt^ from the 
known fluctuation properties of the ideal Bose gas. It is 
also possible, and very instructive, to compute the cy- 
cle probabilities tt^ directly via the infinite-density limit 
of a finite Bose gas (see Section V). Furthermore, we 
will exploit some of the subtleties of the concept of cy- 
cle probabilities to compute the probability distribution 
of the maximum-length cycle at zero temperature, where 
the problem reduces to a study of random permutations 
[5] . Our preceding analysis of the tt^ allows us to under- 
stand why this distribution is essentially unchanged at 
finite temperature, below the condensation temperature 
(Section VI). Finally, we briefly review some known re- 
lations between Bose condensation and the presence of 
infinite cycles for interacting Bose gases (Section VII). 
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II. RECURSION RELATIONS FOR THE 
PARTITION FUNCTION AND THE DENSITY 
MATRIX 

The Landsberg recursion relation [6] gives Zat as a sum 
of only N terms 



1 ^ 



N-kZk, 



(5) 



fe=i 



where the Zk = Zi{kP) = exp {—kl3ea) are the single- 
particle partition functions at inverse temperature k(i. In 
fact, the terms appearing in eqn (5) are the cycle proba- 
bilities: 



TTfe = ZN-kZk/ Zn)- 



(6) 
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FIG. 1: Cycle probabilities TTk and cut-off function for a 
periodic box of size below Tc. 



This is because k particles on a permutation cycle of 
length k at inverse temperature (3 have the same statisti- 
cal weight Zk as a single particle at inverse temperature 
and, furthermore, in the ideal Bose gas, the differ- 
ent cycles of a permutation are statistically independent. 
Thus, in eqn (6) the cycle contributes Zk and the remain- 
ing N — k particles contribute Zjsr-k (see [7]). 
The off-diagonal single-particle density matrix 

PN{r,r',l3)= J da;i . . .dxjv-i 

X p^°''{{xi, . . . ,XN-i,r}, {xi,. . . ,XN-i,r'}, P) , 

is also a sum of permutation-dependent terms, and the 
Landsberg recursion relation can be generalized to a non- 
diagonal single-particle density matrix p(r, r', k(3) rather 
than to the diagonal one: 

PNir, r', /?) = T7 E ^iv-fcP(r, r', k/3), (7) 



fc=i 



^PN{r,r',(3) = Y,^kRT\r,r'), 



N 



(8) 



fc=i 



where V is the volume of the system. In eqn (8), the 
cycle probabilities TTfe are modified by the cut-off function 

i?^"'(r, r'), which is proportional to the ratio of the off- 
diagonal and the diagonal density matrices. In the three- 
dimensional homogeneous Bose gas in a periodic cubic 
box of size I?, we have i?™*(r, r') = p{r,r' ,kl3) / Zk = 
p{0,r-r',kP)/p{0,0,kP). 

As is well known, the condensate fraction is propor- 
tional to the off-diagonal single-particle density matrix 
in the limit |r — r'| oo, a case that corresponds to 
|r — r'l oc L in a cubic box of length L for L ^ oo. The 
cut-off function then vanishes for small k, and terms with 
k^ j0 dominate the sum in eqn (8) [8]. 



III. OFF-DIAGONAL LONG-RANGE ORDER 
AND CYCLE PROBABILITIES 



In this section, we discuss a relation between the off- 
diagonal long-range order parameter and the sum of cycle 

probabilities for the homogeneous Bose gas, with unit 
mass and particle density (m = H = 1, = L^,Tc = 
27r/C(3/2)^/^). To analyze eqn (8), we notice that Zk 1 
for fc > || and Zk - / {2iT(3kf/'^ for k ^. It 



follows from Z]\[-k < Zn that 



< 



l/(27r/3fc) 
1/N 



3/2 



for fc < 1^ 
for > 1^. 



(9) 



We now study the cut-off function i?J:"'(r, r') for |r — 
r'l (X L. For concreteness, wc suppose that the vector 
r — r' has the same components in all three space direc- 
tions: r — r' = L{S, 5, S) with < 5 < ^ and denote the 
corresponding cut-off function and the off-diagonal den- 
sity matrix by i?^"*((5) and by Pn{S), respectively. We 
have 



T 

For k ^ both the numerator and the denominator 

1. For fc < tt^c 



are ~ ^S^/L, so that Rl^'^S) 

numerator is ~ exp {— 3(5^L^/(2/3fc)}, and the denomi- 
nator is ~ 1. This shows that the monotonic function 
is exponentially small for fc <C and equal to 

unity for k (see Fig. 1). The estimates for i?™' 

and for tt^ imply 



hm J2^kRT\5) = 0, 

lim ^ nk[l-RT\S)] =0. 



(10) 
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(see the appendix for a short, but rigorous derivation). 
It follows that 



lim 



N 



/c=constL^ 



0. 



(11) 



This equation relates the probabilities of infinite cy- 
cles to off-diagonal long-range order and (because 
lim|r_r-'|^oo -^PN{r, r', /?) gives the condensate fraction) 
also to the condensate fraction. The relation between 
cycle lengths and off-diagonal long-range order is nat- 
ural because, in order for a cycle to contribute to the 
off-diagonal density matrix for |r — r'| oc L, its de Broglie 
wavelength must be at least comparable to the distance 
i, which means that the particle must belong to a cycle 
of length k with L < y/2Trk(3 (see also [2, 9] for related 
derivations). 



IV. CONDENSATE FRACTION AS A 
DERIVATIVE OF THE CYCLE PROBABILITIES 

In Section III, we studied the sum of the cycle probabil- 
ities. We now discuss the fact that the discrete derivative 
of the cycle probabilities exactly gives the condensation 
probabilities. This relation allows to compute the tt^, 
and it involves neither a cut-off function nor the ther- 
modynamic limit and holds for arbitrary systems of ideal 
bosons. The relation was mentioned very briefly in the 
context of Monte Carlo calculations [4], but it was not 
analyzed in any detail. 

Let us define SN{k, cr) as the canonic partition function 
with N bosons of which exactly k are in state a (see [6]) 
and, similarly, the restricted partition function with at 
least k bosons in the state a, as 



N 



YM{k,a) = J2 SN{k',a) 




■'N-k- 



(12) 



cycle length k 

FIG. 2: Cycle probabilities Hk and their negative discrete 
derivative (not to scale) for = 125 ideal bosons in a periodic 
box of size = 5^ at temperature T/Tc = 0.5 (Af^at = 48.7). 



condensation probabilities to the derivative of the cycle 
probabilities. 

To illustrate the relation between cycles and conden- 
sation probabilities, we consider the example of Fig. 2 of 
N particles in a periodic box, with the cycle probabil- 
ities explicitly computed from the Landsberg recursion 
relation. For k > N/2, the derivative of the cycle proba- 
bilities equals the sum of the probabilities of condensing 
k bosons into state a. More generally, for fc the 
condensation probabilities into the excited states van- 
ish and TTfc — TTfc+i ~ SN{k,0)/{NZj^) directly yields 
the probability distribution of having k particles in the 
groundstate, a distribution with mean value k = (Nq) 
and standard deviation oc L^. On the other hand, for 
small k, as already discussed, the discrete derivative be- 
haves as k~^/'^. It describes the probabilities of condens- 
ing k bosons into excited states. Finally, for <C fc 

and k < (Nq) - L'^/{2T:l3 f/^, the probability of having k 
particles in the condensate vanishes, so that the tt^ are 
constant. 



(5jv(fc, f)/-^Ar is the probability of having exactly k 
bosons in state a.) Summing eqn (12) over all states cr, 
and dividing by N Zm, we recover the cycle probabilities, 



N 



J2YNik,a) = 



-Pkea 



NZ 



-ZN-k = TTfc, 



N 



whose negative discrete derivative, 

= J2SN{k,a)/iNZM), (13) 

a 

thus yields the sum over all a of the probabilities of con- 
densing k bosons into state a. We note that the l.h.s. 

of eqn (13) contains only quantities related to the path- 
integral picture, whereas the r.h.s. involves only energy 
levels. Furthermore, we note that the probability of con- 
densing k particles into excited levels is zero for large 
k, so that eqn (13) effectively relates the groundstate 



V. DIRECT CALCULATION OF CYCLE 
PROBABILITIES 

We have seen in Section IV that the cycle probabilities 
TTfc follow from the known fluctuations of the condensate 
in the canonic ensemble. It is very instructive to com- 
pute them directly via the infinite-density limit for the 
partition function of an arbitrary system. For k ^ j^, 
we have Zk ~ 1, so that TTfc ^ Zn-^/Zn- For illustration, 
the Fig. 3 shows the Zk in a fixed physical system (for the 
homogeneous Bose gas, in a fixed box of size L^), as it 
can be computed numerically from the Landsberg recur- 
sion relation. The limiting value Z^o (at infinite density) 
and the value k = Ng^t of largest variation of Zk can be 
computed exactly. 

Indeed, as the groundstate energy is zero, Zn differs 
from Zn-i only by configurations without any particle 
in the groundstate: Z^r — Z^-i = S'Ar(0,0) (see above 
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FIG. 3: Canonic partition function as a function of A'^ in 
a finite system. 



eqn (12)). Therefore, Zn = EitLi 'S'fe(0,0) with 



Sk{0,0) = 



HE 

.(7>0 7l„-=0. 



exp -p^n^e^ ^E.>o"<..fe- 



cr>0 



Summing this equation over all integers k yields, for ar- 
bitrary states {co-}, 



Zoo — 



J-J- 1 -expf- 



cr>0 



exp (-/Jcct) 



exp(^ 



k 



(14) 



This exact formula for a finite canonical Bose system can 
be obtained without the usual saddle-point integration. 
It agrees with the grand-canonical partition function for 
the excited states at zero chemical potential because at 
high density the condensate serves as a reservoir for the 
excited bosons. 

To determine the value of N for which the partition 
function Zk passes from Zk ^ to Zk ^ Z^ (see Fig. 3), 
we again take a discrete derivative and consider the prob- 
ability distribution of the number of excited particles 

7roo(iVexc = fc) = (Zfe — Zk-l)/Zao, 

whose mean value, the saturation number, is 



-/Vsat = (-^exc) = ^ kWooiNe^c = k) = 



fe=l 



exp {-Peg) 
^^^_-exp(-/3e„, 



= ^fe-l). (15) 



The variance of the number of excited particles is again 
given exactly by the corresponding grand-canonical ex- 
pressions 



1 



^^^exp(/3e,)-l [exp(/?e,)-l]^ 

OO 

= iVsat + ^(fc-l)(zft-l). (16) 



fe=2 



For completeness, we give the expressions correspond- 
ing to eqn (14)-eqn (16) in the homogeneous Bose gas, 
in the limit L ^ oo: Zoo ~ exp [L^/{2nP)-^^^Ci5/2)] 



and, furthermore, N^at/L^ ^ C(3/2)/(27r/3)3/2 ^ .^j^j^ 
r]c the critical density. In addition, the variance of A'oxc 
behaves as _L^/(27r/?)'^ [10], so that the jump of Zk/ZN 
from to 1, for T < Tc, takes place in a window of width 
L^, as already noted in Section IV. 

Putting together the analysis of Zk for a fixed size of 
the system and a varying particle number, and our knowl- 
edge of Zk, we have a complete description of the cycle 
probabilities tt^- Once wo know exactly the behavior of 
the TTft in the thermodynamic limit, we observe that the 
sum of the nk for A; <C ^ gives exactly the critical den- 

sity r]c, and the sum of the tt^ for fc » gives the 
condensate fraction. This conclusion is the same as in 
Section III, but it stems from a microscopic analysis of 
the cycle probabilities. 



VI. MAXIMUM CYCLE LENGTHS 

At finite temperature, the cycle probabilities TTfe are 
easily computed from the Landsberg recursion relation, 
and at T = 0, they are given by TTfc = 1/iV. In Monte 
Carlo calculations the entire permutation of N elements 
can be sampled, starting from the cycle containing the 
element A'', which has length k with probability TTfc. This 
leaves one with a system of N — k particles, for which 
the probabilities of the cycle containing element (N — k) 
can again be computed, etc. (see [7]). (At zero temper- 
ature, the distribution remains Wk = 1/{N — k), at finite 
temperature, it is given by eqn (6).) 

In this section, we use the cycle probabilities ttj. to 
obtain analytical results for p^^^, the probability that 
the longest cycle in a permutation has length k. This 
probability distribution is closely related to the distribu- 
tion TTk both at finite temperature and at T = 0, where 
the ideal Bose gas is equivalent to the problem of ran- 
dom permutations, which has been much studied in the 
mathematics literature. We recover some classic results 
which we generalize to finite temperatures. 

We first note that the cycle probabilities are related to 
the mean number of cycles of length fc, ivk = NiTk/k. For 
k > N/2, TTk coincides with the probability pk that the 
permutation P has at least one cycle of length k. This 
is simply because a permutation of N elements can have 
no more than one such cycle, in other words, because 



^k = Yl 



m=l 



prob. to have 1 _^ 



m cycles of length fc| 



and 



m=l 



prob. to have 
m cycles of length k 



(18) 



and only terms with m = 1 contribute to the above ex- 
pressions if A: > A^/2. 
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Furthermore, any cycle of length k > N/2 must be the 
longest cycle, so that we arrive for all temperatures at: 



■.p^.=nk = Nwk/k (ov k> N/2,yT. (19) 



We next consider the probability distributions pk and 
p^'^'^ for k < N/2 at T = 0, for iV ^ oo. We take the 
continuum limit pk — > p{x) by setting pk = j^p{x = j^) 
(and similarly for p'^^^ix) and tt{x)). 

It follows from the recursive procedure, and from the 
fact that the nk are of order that the probability to 
have more than one cycle of the same length is 0{1/N'^). 
In the hmit N ^ oo, we see that p{x) = fc^x) = 1/x. 

For X < 1/2, the probability p'^'^^{x) is the product of 
the probability to sample x, and the probability that x 
remains the longest cycle in the remaining partition of 
length 1 — a;, in other words: 



prob. that longest cycle in 
[0, 1 — x] is > a; 



At r = 0, and in the limit ^ oo, the probability 
that in a permutation of length A^(l — x) the longest 
cycle's length exceeds Nx equals the probability that in 
a permutation of N elements the longest cycle exceeds 
x/{l — x), and we arrive at 



/" 

Jx/n- 



1- / dx' p'^'^'^ix') 

:/(l-x) 



(20) 



For all g]0, 1[, x < x/{l — x), so eqn (20) can be used 
to compute p'^'^^{x) from x = 1 downwards to arbitrary 
precision (see Fig. 4). Alternatively, we can transform 
eqn (20) into a differential equation 



X 



^(x) 



X(l-X)2^ 



1 -X 



Explicit formulas for p™^'^(x) are easily obtained, starting 
by inputting, from eqn (19), p™^''(x) for i < x < 1 to 
obtain p"^^ in the window x e [|, etc, 



-{x) = 



for 



< X 



^i[l-log(i^)]for 



3 < X < 2 , 



(21) 



but they become cumbersome (see also [11]). 

At finite temperature, the cycle of the element N is 
sampled from a distribution tt^, as shown in Fig. 2 rather 
than from the flat distribution. However, the longest 
cycle is again of length oc A'', and it is sampled from 
an essentially constant distribution in the interval x € 
[0, (A^o) /N]. In the thermodynamic hmit, for T < Tc, 
the distribution p™'^^{x) of the properly rescaled variable 
X = k™^^ / (No) trivially agrees with the one obtained at 
T = 0. In particular, the mean length of the longest cycle 
is given by 0.624 (A^o) at all temperatures T < [5]. At 
sufficiently low temperature, the rounding at x ~ 1 is 
exactly described by eqn (19) (see Fig. 4). 



T=TJ2 iV-512 
n A=512 
T=0 



a, 

>5 



03 

o 




cycle length x = ^"''^/(Aq > 



FIG. 4: Probability p'^^^{x) from a numerical integration of 
eqn (20) at T = 0, compared to the distribution p™^'' obtained 
from a Monte Carlo sampling at T = Tc/2, and the exact 
expression of eqn (19) for large a;, at T = Tc/2. 



VII. INTERACTING BOSE GAS 

For the interacting Bosc gas, the relation between off- 
diagonal long-range order and the infinite cycles has been 
studied for generalized mean- field models ([2, 12]), but 
it has proven difficult to establish rigorous results for 
models with realistic interactions. In the general case, 
the definition of the superfluid density in terms of the 
winding numbers [13] proves that for the homogeneous 
Bose gas in three dimensions, the presence of infinite cy- 
cles is equivalent to a non-zero superfluid density. It is 
generally admitted that in homogeneous systems, Bose- 
Einstein condensation implies superfluidity, so that the 
presence of a condensate would imply the existence of 
infinite cycles in those systems. 

On the other hand, the direct link between the cycle 
probabilities and the condensate fraction, as described in 
Section IV, cannot hold for interacting systems. This is 
easily seen at T = because the groundstate of any sys- 
tem of identical but distinguishable particles has bosonic 
symmetry. Therefore, the distinguishable-particle den- 
sity matrix 



Pi{xi 



,xn},{xp^,. .. ,xp^-},l3 = oo) 



is independent of the permutation P. It follows that at 
T = 0, and for finite N, all permutations have the same 
weight and the cycle probabilities again satisfy tt^ = 
At the same time, the condensate fraction of interacting 
systems at T = is less than one[14], and the distribution 
of the TTk cannot be related to the distribution of Nq in 
the same manner as in the ideal Bose gas. 



VIII. CONCLUSION 

In conclusion, we have discussed the relation between 
the off-diagonal long-range order, the condensate frac- 
tion, and the cycle probabilities in the ideal Bose gas. 
Only long cycles (fc > constL^ in the homogeneous case) 
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contribute to the off-diagonal long-range order parame- 
ter, and Bose Einstein condensation is equivalent to the 
presence of infinite cycles. We have also discussed the 
probability distribution for long cycles, and the general, 
non-intuitive, link between the integer derivative of the 
cycle probabilities and the number of condensed bosons. 
This integer derivative provides us with a purely topo- 
logical characterization of the condensate fraction in the 
ideal Bose gas. Our knowledge of the cycle probabilities 
allows us to study the distribution of the longest cycle. 
We have remarked that the <T < case can be un- 
derstand from the T = case. The mean length of the 
longest cycle is proportional to the number of condensed 
bosons. For the interacting Bose gas, at the present time, 
only the winding-number formula[13] has been rigorously 
shown to relate the topology of Feynman paths to an- 
other characteristic of interacting Bose-condensed sys- 
tems, namely the superfluid density. The link between 
Bose-Einstein condensation and superfluidity will need 
to be explored further in order to better understand the 
cycles of the interacting Bose gas. 
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APPENDIX A: DERIVATION OF EQN (10) 

In this appendix, we prove eqn (10). We consider the 

sums separately in the regions a) d) indicated in Fig. 1 
and show that they all vanish individually. 

In region a), 7rfci?^"*((5) < tt/c exp [— 3(5^LY/?] , and 
the normalization of the TTfc shows that the partial sum 
vanishes in the limit L ^ oo. In b), TTfei?^"' < 
7ri2-ei?^"' < const/L^'^^"'^)/^ and the partial sum is 
smaller than constL^^/^"^. In c), as 1 - R%^^{5) de- 
creases with k, TTfc [1 - i?^"*((5)] < TTL2 [1 - R'j^^iS)] and 
the partial sum is dominated by constL*^"^. Finally, in 
d), TTk [1 - -RrH'^)] < TTfc [1 - i?^f+.((5)] . The normal- 
ization condition of the tt^ and the relation ^ 
imply that this partial sum vanishes for L ^ oo. With 
a suitable choice of e (for example e = |), these sums all 
vanish for L — ^ oo. In fact, any function e(L) ^ 1/logL 
that satisfies e(L) < 2/3 for L > Lq can be used. 
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